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ABSTRACT 

Context. Magnetic clouds (MCs) are "magnetized plasma clouds" moving in the solar wind. MCs 
transport magnetic flux and helicity away from the Sun. These structures are not stationary but 
feature temporal evolution. Commonly, simplified MC models are considered. 
Aims. The goal of the present study is to investigate the dynamics of more general, radially 
expanding MCs. They are considered as cylindrically symmetric magnetic structures with low 
plasma f}. 

Methods. The self-similar approach method and a numerical approach are used. 
Results. It is shown that the forces are balanced in the considered self-similarly evolving, cylin- 
drically symmetric magnetic structures. Explicit analytical expressions for magnetic field, plasma 
velocity, density and pressure within MCs are derived. These solutions are characterized by con- 
served values of magnetic flux and helicity. We also investigate the dynamics of self-similarly 
evolving MCs by means of the numerical code "Graale". In addition, their expansion in a medium 
with higher density and higher plasma f3 is studied. It is shown that the physical parameters of the 
MCs maintain their self-similar character throughout their evolution. 

Conclusions. A comparison of the different self-similar and numerical solutions allows us to con- 
clude that the evolving MCs are quite adequately described by our self-similar solutions - they 
retain their self-similar, coherent nature for quite a long time and over large distances from the 
Sun. 

Key words. Magnetohydrodynamics (MHD) — Magnetic fields — Plasmas — Sun: solar 
wind 

1. Introduction 



It is well-known that coronal mass ejections (CMEs) are one of the most significant forms of 
solar activity. They carry enormous masses of plasma threaded by the magnetic field away into the 



2 Dalakishvili et al.: Investigation of Dynamics of Self-Similarly Evolving Magnetic Clouds 

interplanetary medium. Further away from the Sun, these large-scale, dynamical plasma structures 
are commonly called interplanetary coronal mass ejections (ICMEs). Magnetic clouds (MCs) form 
a subset of ICMEs (Klein & Burlaga [T982l Burlaga fl99T1 Farrugia et al. 119951 . Spacecrafts crossing 
the central parts of such MCs provide valuable information about their physical characteristics. It 
turns out that MCs have a strong magnetic field, low proton temperatures (low plasma f3, compared 
to the ambient solar wind with the same speed) and they feature a substantial and smooth rotation of 
the magnetic field vector. These three features of MCs are selected as signatures of MCs (Nakwacki 
et al. 120081 ) . The MCs are also characterized by a coherence of the magnetic field (low level of 
fluctuations). The radial dimension of a MC is typically « 0.25 AU (at 1 AU). 

These in situ observations of the physical properties of MCs are considered as important steps 
towards the prediction of the geophysical effectiveness of their interaction with the Earth's magne- 
tosphere, for space weather forecasts and related issues. 

Different models for the structures of magnetic clouds have been proposed. There is no gen- 
eral agreement about the large scale structure of MCs. Commonly, the local structure of MCs is 
considered in the form of cylindrically symmetric force-free configurations (Burlaga 1988 [19911 
Demoulin & Dasso 2009). It is often suggested that the ends of MCs connect to the surface of 
the Sun while, according to other models, MCs are described as tori (Vandas et al. 2006, 2009, 
Romashets et al. 120061 120071 ). In a number of studies, MCs are considered as force-free, static, axi- 
ally symmetric flux ropes and their magnetic field is constructed on the basis of Lundquist's model 
(Burlaga 1988 Lepping et al. 119901 Farrugia et al. 119931) . Observations show, however, that MCs 
do not stay static but expand while propagating in the solar wind and they keep expanding well be- 
yond 1 AU (Burlaga fTWTl Demoulin[2008] Demoulin & Dasso l2009l Bothmer & Schwenn [l998b . 
In a large majority of the cases it is observed that the frontal parts of the MCs propagate with higher 
velocities than their back regions. This shows that, with respect to the MCs own cylindrical set of 
coordinates, the radial size of those cylindrical MCs increases (Nakwacki et al. 120081 . Theoretical 
models including the effect of radial expansion have been proposed before (Osherovich et al. 1 19951 
Farrugia et al. 1993 Nakwacki et al. 120081 ). In these models, only the radial expansion is taken into 
account and solutions have been found for all plasma parameters. There are other studies (Shimazu 
& Vandas 2002 Demoulin & Dasso 2009), however, where the axial expansion is also included. 

Previous studies also showed that inside MCs the density drops as dr 1A (Bothmer & Schwenn 
119981 . i.e. the volume of MC increases as d 2A , where d denotes the distance from the Sun. The 
radius of the MCs, denoted by (R), also increased and at a rate changing with the distance, viz. as: 
R ~ d * (Bothmer & Schwenn |1998l ). Since the surface of the MCs cross section perpendicular to 
its axis increases as R 2 ~ d 16 , and the MCs volume increases as d 2A , the MCs longitudinal size 
should increase as d & . Therefore, according to Bothmer & Schwenn's data, the MCs are radially 
expanding and also show an extension along their axis. 

In the present study, we consider self-similarly expanding cylindrical MCs that are able to 
expand both in the radial and longitudinal directions. We consider the problem in the frame of 
the MC and in cylindrical coordinates related with the MC, i.e. with a longitudinal axis Z that 
coincides with the MCs axis. Overall cylindrical symmetry of the MC is assumed. Based on these 
assumptions, we derive the appropriate full set of non-stationary MHD equations and find their 
analytical solutions. The logical and natural consequence of the assumptions of self-similarity and 
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cylindrical symmetry is that the dynamic forces acting upon the MCs are balanced. The solutions 
include expressions for the plasma magnetic field, velocity, mass density and thermal pressure. 

An important feature of our model is that certain significant characteristics of the MCs — 
magnetic flux and helicity — are conserved. We separately consider also the particular case of a 
MC that is allowed to expand only in the radial direction. It can be shown that in this case, the 
MHD equations do not have any physical, self-similar solution. 

2. Self-similar expanding MC models 

2. 1 . General equations and self similar expansion 

In order to perform an analytic study of the dynamics of magnetic clouds, we have to start from the 
full set of MHD equations: 

V-B = 0, (1) 
d,B = V x [V x B] , (2) 
d,g + V • (g\) = 0, (3) 
g[d, + (V • V)]V = (1/4tt)(V x B) x B - V • p, (4) 

In these equations, p denotes the thermal plasma pressure, g is the density, V is the velocity 
field and B denotes the magnetic field. 

In a number of previous studies, the MCs were considered as cylindrical magnetic structures, 
characterized by axial symmetry. In the present consideration, both symmetry along the Z axis 
{d, = 0) and the azimuthal symmetry (d^ = 0) are assumed. The axially-symmetric magnetic field 
can then be expressed in the following way 

B = [0, B v , BJ, (5) 

where B^ = B^(r, f) and B, = B z (r, f). Note that this representation satisfies the solenoidal condition. 

The self-similar approach, adopted here, implies that the temporal evolution of the physical 
functions is controlled by the following self-similar variable: 

where €>(f) denotes a function of time. Let us search solutions of the MHD equations in the follow- 
ing form (in analogy with Low |1982j l: 

By = ® S Qv(&, (7a) 
B z = ^QMl Ob) 

g = <s>°p({), (7c) 

p = cw) 

One can see that, the type of solutions introduced by Eqs. (7a-7d) evolve self-similarly and are 
characterized by a particular time scaling. 

Here Q v , Q z , p and p are functions of the self similar variable <S) S , O ", cD a and show 
the time scaling of the azimuthal and longitudinal components of the magnetic field, the plasma 
density and the plasma pressure, respectively. 
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2.2. Solution of the induction equation 

We consider both a radial and a longitudinal expansion of the MC but no motion in the azimuthal 
direction is considered. In this case the Eulerian velocity field of the plasma, V, can be expressed 
in the following way: 

V = [V,, 0, VJ. (8) 

Here, we assume that the radial component of the velocity V r = V r (r, t), and the z-component 
V z = V z (z, t), i.e. we assume that the MC maintains its cylindrical shape during its evolution. 
After substitution of Eq. (8) in Eq. (2) we derive: 

d,B- + -d r (rV r B~) = 0, (9a) 
r 

After taking into account relations (6), (7b) and the relations (Al) and (A2) given in the appendix, 
Eq. (9a) can be rewritten as follows: 



Qz 



V r 

O-O + — + (bd r V r 



+ Q' z [v r ~Z<b] = 0. (9b) 



Here Q' corresponds to dQ z (^)/d^. Equation (9a) (therefore equation (9b)) is satisfied for arbi- 
trary Q- only when: 

V r -£® = 0, (10a) 

and 

o-O + -j + Od r Vr = 0. (10b) 

From Eq. (10a) and Eq. (10b) follows that radial component of the Eulerian plasma velocity is 
described as follows: 

V r = r®/Q>, (11) 

and 

o-=-2, (12) 

where <& is the function of time mentioned in Eq. (6). 

One can check that for <j - -2 the longitudinal magnetic flux <f>. is conserved. Nakwacki et al. 
(2008) analyzed different MC models and derived expressions for the magnetic flux, the magnetic 
helicity and the magnetic energy per unit length along the flux tube. The models which are in good 
agreement with observations are characterized by the conservation of <p z , see also Berdichevsky et 
al. (120031 

Let us analyze the ^-component of the induction equation, Eq. (2): 

d t (B v ) + B^Vz + dJYrBJ = 0. (13) 
The combination of Eqs. (11) and (13) leads to the following important relation: 

(6 + 1)6/0 + d z V z = 0. (14) 



After taking into account expressions (5) and (8) in combination with the assumption of az- 
imuthal symmetry, one can see that radial component of the induction equation, Eq. (2), is auto- 
matically satisfied and does not lead to any additional restrictions. 
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2.3. Self-similar solutions 



After inserting the expression for the plasma density (7c), together with the velocity from Eq. (8), 
with Eq. (11) for the radial component, in the mass conservation law Eq. (3), we obtain another 
important relation, viz. 

(or + 2)0/«D + 8 Z V Z = 0. (15) 

Obviously, in order to have consistency between Eqs. (14) and (15), one should have: a + 2 - 
6+1. 

The z-component of the equation of motion, Eq. (4), helps to derive an expression for the z- 
component of the plasma velocity: 

d t V z + V z d z Vz = 0. (16a) 

Let us try to solve the partial differential equation (16a) by using the variable separation tech- 
nique, i.e. we assume that 

V z (z,t)=Z(z)T(t). (16b) 
Substitution of expression (16b) in Eq. (16a) yields: 

fZ + Zr 2 Z'=0, (17a) 

here F = d,F denotes the first order time derivative of a function F. Hereafter we will use, for 
indicating second order derivatives, the notation: F = djF. While Z' stands for dZ/dz. 
It follows from Eq. (17a) that: 

f 

-— = Z' = const, (lib) 
Equation (17b) can be decomposed into two ODEs, viz. 

-|r = *> (17c) 

and 

Z' = A. (lid) 

Here, A is an arbitrary constant. 

After solving the ODEs (17c-17d) with the assumption that, at the surface z = 0, V z = 0, we 
derive the following expressions: 

T=- ? — , (18a) 

l+AT t 

and: 

Z = Az, (18i) 

where A and Tq are constants. 

After inserting Eqs. (18a-18b) in Eq. (16b), we obtain the wanted expression for V z : 

v - = TTi? <1!> > 

Here, k = ATq. 

We assumed that locally the MC could be described as a cylindrical structure. Let us investigate 
the evolution of the length L of this cylindrical structure. For this purpose, let us describe the 
temporal evolution of the z-coordinate of the plasma element located at the position z = L at time 
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f. The Lagrangian velocity of this element coincides with the Eulerian velocity of the plasma flow 
at time f and z — L. If at a certain time the coordinate of this element is L, then its Lagrangian 
velocity is: 

dL 

V L =—. (20a) 
at 

From Eq. (19) we then have: 

f = T~T' ™ 
dt l+kt 

The solution of this ordinary differential equation (20b) gives the following expression for the 
longitudinal size of the considered cylindrical structure: 

L = Lo(l+kt), (20c) 

where Lq is the length of the cylinder at t = 0. 

The radial component of the equation of motion, in combination with the expressions for the 
magnetic field, the velocity and the plasma density leads to: 

<£^i>p = F r . (21) 

Here, F r denotes the radial component of the total force. In terms of p, Q VtZ , and this force 
component can be expressed as (details of the derivation are given in the appendix): 



1 

F r = 

4tt 



1 (O 2 

0)5 



S^-y , (22) 



where QL Z - dQ^-Jdt;, and p' = dp/di;. In order to have a self-consistent time scaling for all terms 
in Eq. (22), one has to require that 6 — -2 and (3 = -4. From a comparison of Eq. (14) to Eq. (15), 
follows that if 5 = -2, then a = -3. 

At the same time, from Eq. (15) and Eq. (19): 

O k 



» l + kf (23) 
Equation (23) is an ordinary differential equation in terms of 0(f). After solving this ODE, we find 

the following explicit expression for 0(f): 

= o (l+/fcf). (24) 

Here 0o is a constant parameter. 

The substitution of expression (24) in Eq. (21) leads to an important conclusion: calculating 
the magnetic and pressure gradient forces, we see that for the self-similarly evolving, cylindrical, 
axially-symmetric structure the magnetic force, F m = 1/(4^)(V x B) x B and the thermal pressure 
gradient force, F p = - V -p are exactly balanced, i.e. 

F r = F,„ + F p =0. (25) 

If we associate the value £o of the self-similar variable £ with the boundary of the MC, then the 

expression of the MC Lagrangian velocity is given by (Low |1982l ): 

dR c/0 

V* = - t =^- t . (26) 
After substitution of the expression (24) for in Eq. (26), we can derive a time-dependent solution 
for the MC radius: 

R=R (l+kf). (27) 
Note that the form of this expression coincides with the one given by Nakwacki et al. (2008). 
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2.4. Plasma and force-free field evolution 



The rest of the solutions readily follows from the derived equations, yielding: 



V r = ^ (28a) 



~ (28fo) 



(1+A:/) 3 
and 

P = — 7. (28c) 

Here, g and p are arbitrary functions of £ = r/O. 

After analysis of the expressions for pressure and density (28b-28c), one can check that, for 
systems characterized by entropy conservation, the entropy conservation law is satisfied only if the 
poly tropic index y = 4/3. Actually, this is a common feature of all different self-similar systems 
(Low [1982l Farrugia et al. [19951 Finn et al. 120051 . 

From various observations it is known that MCs are characterized with low plasma yS's (Burlaga 
et al. 1981 Burlaga 1991 Bothmer & Schwenn 1998). The thermal pressure term in the total force 
could be neglected and this implies that we have to construct a force-free magnetic field that evolves 
in a self-similar way. The cylindrically symmetric force-free structure of the MCs magnetic field 
is indeed advocated by a number of researches (Burlaga 1988 Lepping et al. 1990, Farrugia et al. 
[19931 Farrugia et al. [19951 Nakwacki et al. 120081 Demoulin & Dasso |2009l ). A force-free magnetic 
field satisfies the following relation: 

V x B = yuB. (29) 

If we rewrite the vectorial equation (29) for each component of vectors, taking in to account ex- 
pressions (7a), (7b) and (12), we obtain: 

-fi' = yfS>Q v , (30a) 

e; + y = ^q z , (30*) 

here Q' , stands for dQ^/dg, If we take the derivative of both terms of Eq. (30a) with respect to 
the variable we get: 

-Q' z ' = f i®Q' ip . (30c) 
Here it was assumed that p. does not depend on In general, however, /j could be a function 

If we take in to account expressions (30a) and (30c), we can derive from Eq. (30b) an ordinary 
differential equation for Q z : 

Q'l + Q-+i/® 2 Q z = 0. (31a) 
With the following transformation of variables: x = n$>%, we can rewrite Eq. (31a) as follows: 

d 2 Q z 1 dQ z 

-T^ + --T 1 +e 2 = 0. (31*) 

dx L x dx 

Actually Eq. (31b) is a Bessel equation of zero order, with the following solution: 

Q z = J (x) = C ./o(/^), (32a) 
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where Jq(x) is the Bessel function of the first kind, Co is a constant parameter. Notice that the 
solution which is not characterized with a singularity at x = has been chosen. The substitution of 
Eq. (32a) in relation (30a) leads to an expression for Q^: 

Op = Ji(x) = C /i (/!*£), (32fo) 

with J\ (x) the Bessel function of the first kind. 

From Eqs. (32a-b) we see that and Q z are the functions of yu<l>£. Since we assumed 
above that /j is not a function of £ that $ is only a function of t (Eq. (6)), and that Q^ z are 
functions of only £ (Eqs. (7a,7b)), it follows that yuO = const. The substitution of Eqs. (32a), 
(32b), (24) and (6) in Eqs. (7a) and (7b) respectively, taking into account that 6 = a = -2, 
leads to the following expressions for the components of the magnetic field: 

B r = 0, (33a) 

Ba, = — ^-^Jii 1, (33b) 

* (l+kt) 2 \r (l+kf))' 

and 

Bz = ,« B ° , /o I ,/ ,J . (33c) 



(l+kt) 2 \r (l +kt) 
where B and r are constants. 

(Co/^y is changed by B Q and /U^/O is substituted by l/r .) 

From Eqs. (33a - c) we can calculate important expressions for the magnetic flux and the 
helicity (Nakwaci et al. 2008) associated with the MC: 



2tt r B () 
X (l+kt) 2 



= — R — hs J i(X«), (34a) 



*(l+£f) 2 

and 

H = —R 2 L(j](xR) - Jo(xR)Ji(xR) + J 2 (*R)), (34c) 
X (l + K») 

where ^ = l/(r (l + &))• 

By its physical meaning O z is the magnetic flux across the surface perpendicular to the axis of 
a MC, while cD^ is the magnetic flux across the surface defined by the magnetic axis and the radial 
direction. Moreover, R denotes the radius of the MC and L is longitudinal length of the cylindrical 
structure. By inserting in Eqs. (34a - c) the corresponding expressions for R and L we find: 

cD, = 27iRoroBoJi { — - | = const, (35a) 



% - B r Lo 

and also 



ro 



= const, (35b) 



H = 27rr R 2 B 2 L 



= const. (35c) 



ro I \ r I \r / \ r 
From these results it follows that the obtained solutions ensure the conservation of magnetic 
flux and helicity inside the cylindrical MC described by our model. 
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3. Radially expanding MCs 

The purpose of this section is to find solutions for the physical variables in the case where only 
the radial size of the MC increases. One can see that the solutions in this case do not remain self- 
similar, although initially a self-similar expansion is assumed in the radial direction. 

Below we consider MCs that are expanding only radially; i.e. with V z = 0. In this case Eq. (14) 
implies <fc = 0or<5=-l. The case with <t = corresponds to the stationary state, which is trivial. 
Let us consider the case when <fc + but 6 = -1. In order to provide a consistent time-scaling of 
all terms in Eq. (22), we have to satisfy: 

e^e; + y = ^fmQJ = o. 0&1) 

An analysis of Eq. (6) = r/<D(f)), Eq. (7a) (fi v = O' 5 £></;(£» and Eq. (36a) leads to the 
following expression: 

-d r {rB v ) = 0. (36*) 
r 

Here, it is taken into account that d r = 3^/0(f). 
From Eq. (36b) we can conclude that: 

B<p = — , (36c) 
r 

with C = const. 

Note that the expression for B v is characterized by a singularity at the axis (r = 0). It seems 
reasonable to conclude that, if we do not consider the axial stretching of self-similarly evolving 
MCs, we can not obtain a physically valid solution for the B v — component on the axis of the MC. 



4. Numerical study: higher density and higher plasma j3 case 

In this section, we investigate the evolution of MCs in a medium by means of the model described 
in Section 2. For this purpose, the Lagrangian numerical MHD code "Graale" (Finn et al. 2004) 
is used, which enables us to check whether the above-obtained solutions maintain their self-similar 
nature when they propagate in a medium. In the numerical code, azimuthal and cylindrical sym- 
metries are implied. Furthermore, it is assumed that the magnetic structure expands uniformly in 
the longitudinal direction, in other words V z = zL/L; with L the length of the cylinder. An analy- 
sis of Eqs. (20a-20c) and Eq. (19) shows that the derived expression for the longitudinal velocity 
coincides with the one implemented in the code. 

The assumptions of cylindrical and azimuthal symmetries in the code and the prescription of 
the character of the longitudinal motions makes the numerical simulations ID. In the numerical 
runs the units of the physical parameters are chosen as follows: the unit length L,„„, = Q.IAU = 15 • 
10 6 km is of the order of the MCs radius at 1 AU, the unit magnetic field B mlt = 3 nT, and the unit 
number density «„,„■, = 10 cm -3 . After taking into account that the proton mass m p « 1 .7 • 1 0~ 27 kg, 
one derives that the unit mass density is p„„„ = m p rio - 0.8 ■ 10~ 14 kg/m 3 . The unit speed is the 
Alfven speed corresponding to p unlt and B mll \ V m it - Voa = 20.5 km/s, which is of the order of 
the MCs edges expansion velocity in the frame of the MC (Vandas et al., 2009), and the unit time 
tunit = jr 1 = 200 h. A domain with R mm = and R max = 10 is discretised with 2000 grid cells. The 
time step used in the simulations is At = 5 • 10~ 7 (/?„„■„, R max and At are given in units introduced 
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above-L„„ if and t u „i,). Open boundary conditions are applied. Inside the calculation domain, we 
introduce initial conditions for the physical variables in two different regions: inside the MC and 
outside the MC. 

Evidently, the solutions inside and outside the magnetic structure should satisfy the following 
jump conditions across the surface of any MC: 



[gv r ] = (37a) 
) (37b) 



Q^r +P+ ^-B 2 



1 3 +( _X^ 

2 y-l w 



= (37c) 



and, finally, 

[gv r v,] = (37a!) 

Here, [■] denotes the jump of the quantity between the brackets across the surface of the MC. 
Also, v r = V r - V s , where V r and V s are the plasma and the MC's surface velocity, respectively, 
while v, denotes the plasma velocity tangential to the surface of the MC. Equation (28a) and 
Eq. (26) show that v r = 0, which is logical for ideal MHD. Equations (37a,37c,37d) are satis- 
fied for arbitrary values of the plasma density and Eq. (37b) leads to the condition: 



1 3 



= 0. (38) 



We know that the plasma mass density inside the MCs is lower than outside them and the 
plasma B within a MC is lower than in the ambient plasma. We therefore consider B ~ 1 in the 
ambient environment and B ~ 0.1 inside the MC. For the magnetic field within the MC, we use 
the solution expressed by Eqs. (33a-33c). For the magnetic field outside the MC, we assume that 
the azimuthal component of the magnetic field B^,, = 0, while the longitudinal component B zou , 
is uniform. We also assume that the mass density and the thermal pressure are uniform in both 
regions of the computational domain. 

Bearing in mind these assumptions and jump conditions (37 - 38), one can find explicit expres- 
sions for the plasma pressure and magnetic field outside the MC. 

Figure 1 represents the numerical solutions for the plasma mass density and velocity, while 
Fig. 2 shows solutions for the magnetic field components at different moments in time. On panels a 
and b we plotted the dimensionless values of the density, velocity and magnetic field. On panels c 
and d the dependence of the modified values of the physical parameters on the self-similar variable 
is presented. 

The dependence of the modified mass density and velocity as well as the dependence of the 
modified magnetic field components on the self-similar variable clearly shows that our solutions 
maintain their self-similarity in the course of the MC expansion. 



5. Discussion and Conclusions 

In this paper, we presented a detailed derivation of a class of self-similar analytic solutions of the 
MHD equations for both radially and axially expanding MCs and a numerical investigation of these 
solutions. The usage of the self-similar approach is quite common for the modeling of various kinds 
of solar plasma structures, flows and eruptions (Low 1982, Osherovich 1993, Farrugia et al. 1 19951 
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Nakwacki et al. 120081 Shapakidze et al. 1201 01) . In most of the previous studies, however, only the 
radial expansion of the MCs was considered. In the present study, we took into account also the 
axial stretching of the MCs, which is a common observed feature of at least some MCs. We have 
obtained explicit analytical expressions for the magnetic field, the plasma velocity, the density and 
the plasma pressure. Essentially, our solutions maintain their self-similar nature during the whole 
course of their evolution and propagation through the solar wind. These solutions are complete and 
well-defined, fully analytic and, moreover, in the particular case of the absence of the longitudinal 
expansion, our solutions self-consistently match the analytic solutions derived by other authors 
(Farrugia et al. [19951 

Note that for the class of solutions introduced by Eqs. (7a-7d), the assumptions of self-similarity 
and axial and azimuthal symmetry lead to the fact that i>(f) is a linear function of time (<t>(f) is the 
time dependent function of the self similar variable £ = r/<D(f)). In this case, the forces within the 
MCs are bound to be balanced. We can thus conclude that the case in which the magnetic structures 
are characterized by a low plasma-/?, corresponds to the force-free magnetic field case. Remark that 
this result is also in agreement with the conclusion of previous studies. We therefore believe that 
this is a correct and proper time-dependent generalization of the widely used stationary Lundquist 
model (Lundquist 1950). Note also that in their recent papers Vandas et al. (2006, 2009) made 
a comparison of the generalized Lundquist model with observations and found good agreements 
between this classic model and the experimental data. 

It must be emphasized that our study is not the only one in which the axial stretching of the 
MCs is taken into account together with their radial expansion. As a matter of fact, Shimazu & 
Vandas (2002) also considered MCs with similar properties. In this particular paper, the authors 
used the mathematical approach introduced by Osherovich et al. (1995). In order to separate the 
time-dependent parts of the solutions from multiplicative functions of the self-similar variable only, 
Shimazu & Vandas (2002 ) imply a so-called "separable magnetic" field, which was introduced in 
Osherovich et al. ([T995). The approach introduced in Osherovich et al. (119951 ). turns out to be 
quite restrictive because it requires an ad hoc relation between the different components of the 
magnetic field. Also, in order to separate the time-dependent part from the coordinate dependent 
parts in the momentum equation, in addition to the polytropic law, the authors introduced a specific 
mathematical expression for the thermal pressure (see Eq. 17 Osherovich et al. 1995 ). Actually, 
the mentioned expression relates pressure and mass density (see Eqs. (13-17) Osherovich et al. 
119951 ). In our study, on the contrary, we used B.C. Low's approach (1982) and required a similar 
time-scaling for all parts of the Lorentz force and the force caused by the gradient of the thermal 
pressure. We argue that our approach is more general and puts less non-physical restrictions upon 
the physical parameters. 

Another difference of the results presented here with those of Shimazu & Vandas (2002 ) is with 
the temporal expansion scaling. Shimazu & Vandas assume that the longitudinal and the radial 
expansion have the same time scaling, while in our work this is not assumed but it rather logically 
follows as the by-product of the accurate solution of the MHD equations. In their paper, the time- 
dependent function of the self-similar variable is characterized by a linear dependence on time 
only when the thermal pressure is zero, while we have derived an explicit expression of this time- 
dependent function of the self-similar variable and it has been shown that <£ is a linear function of 
the time variable and it does not depend on the character of the pressure function. 
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Yet another difference between the results presented here and the results of Shimazu & Vandas 
is related to the structure of the MC magnetic field. Shimazu & Vandas, in order to derive explicit 
expressions for the magnetic field, used the assumption that the magnetic structure of the MC 
is described by a force-free magnetic flux rope. In our study, however, we derived an explicit 
expression for the magnetic field. We have derived ordinary differential equations (Eqs. (30a)- 
(30b)) for the functions describing the components of the magnetic field, after solving the equation 
of motion for the case which corresponds to a low plasma (3 within the MC. For a particular type 
of parameters, we have found explicit, analytical solutions for the components of the magnetic 
field (Eqs. (33a)-(33c)). Note that these expressions are a time dependent generalization of the 
well-known Lundquist solutions (Lundquist 1950, Burlaga 1988). 

Our model implies the conservation of magnetic flux and helicity by design, which is satisfac- 
tory and in good agreement with previous investigations (Nakwacki et al. 2008 Demoulin & Dasso 
120091 Kumar & Rust fl996t . 

For further confirmation of the validity of our solutions, we investigated the dynamics of mag- 
netic clouds numerically. In the numerical code Graale we introduced our self-similar solutions as 
initial conditions. The obtained numerical results showed that during the evolution and propagation 
of these MCs, their physical variables maintained their self-similar character. This circumstance 
was illustrated by Fig. 1 and Fig. 2. 

Obviously, the class of solutions found in this paper is quite idealized. The assumptions about 
the self-similar evolution and the consideration of a cylindrical symmetric structure are quite well- 
justified, but real MCs show self-similar coherence and cylindrical symmetry only approximately. 
Hence, in a future study it would be reasonable and interesting to consider more realistic config- 
urations. There are several issues related to the model which can be tested and generalized in a 
forthcoming study: 

1. Our assumptions, just like in previous investigations (Low 1982, Farrugia et al. 1 19951 Finn et 
al. 120041 Shapakidze et al. 1201 Oi l, for the systems where entropy is conserved, put a restriction 
on the value of the poly tropic index y = 4/3. We would like to develop a model that helps to 
avoid this restriction. 

2. Our model describes the plasma dynamics only inside the MC. In the near future, we plan to 
investigate the interaction of an MC with its environment by constructing consistent solutions 
of the MHD equations outside the MC. 

3. We investigated the obtained analytical solutions numerically using a ID MHD code and a sim- 
ple model for the flow outside the MC was implemented. It would be interesting and reasonable 
to study the MC evolution also with 3D numerical codes, where a more complicated and re- 
alistic background flow can be implemented (in preparation). For this purpose the obtained 
solutions could be used as initial state in the 3D numerical simulation codes. We are interested 
in an investigation of the different possible boundary conditions on the surface of the magnetic 
cloud. 
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Appendix A: 

In this appendix we would like to give some details of derivation for Eq. (21) and Eq. (22). 
During the process of derivation were implied following relations: 

d r = ±d € , (Al) 

d, = -|*fl f . (A2) 

These relations follow from Eq. (6). 

Equation (21) is the radial component of the equation of motion-Eq. (4). In cylindrical coordi- 
nates the left term of Eq. (4) can be written as follows: 

p [8, + (V • V)V r ] = p [8,V r + V r d r V r ] . (A3) 

If we take into account relation (6) and expression for V r (Eq. 11) we can obtain following 
expressions: 



* V, ' = r *- r i*j' (A4) 
V r d r V r = r(^\ . (A5) 



After substitution of Eq. (A4) and Eq. (A5) into Eq. (A3) we get: 

p[3, + (V-V)V r ]=p£*. (A6) 
If we combine of Eq.(7c) and Eq (A6) we obtain: 

p[d t + (y-V)Vr\ = <V&p, (A7) 

the left term of Eq. (21). 

In order to derive the first part of the right term of Eq. (22) let us introduce the following 
notation: 

VxBiJ. (A8) 
After taking into account expressions for magnetic field (7a-7b) with <x = -2 we get: 

J r = X -d^B z - 8^ = 0, (A9) 
J 9 = d z B r - d r B z = -d r B z = -||, (A10) 

and 



.S-l 



J z = -d r (rB v ) - -d v B r = -d r {rB v ) = ^— + Q'fl s l . (All) 
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Since we know expressions of vector J we can derive the vectorial product of J and B which 
represents the first part of the right term of Eq. (22). 

Combination of relation (Al) with the expression for pressure Eq. (7d) leads to the expression 
of the second part of the right term of Eq. (22). 
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Fig. A.l. Snapshots of the evolution of the plasma density and velocity field. Panel a) represents the 
dependence of the plasma density on the radial coordinate for four moments in time; Panel b) shows 
the dependence on the plasma velocity on the radial coordinate for four time moments; Panel c) and 
Panel d) illustrate the dependence of the modified density and the modified velocity, respectively, 
on the self similar variable. Parameter values for this case are k — 2, Bo = 1, Qout = ^>Qm, Pin =0.1, 
p out = 1, and (Do = 1. Black line corresponds to t — 0, red line represents the moment t = 0.5/3, 
green line shows the time moment 1 /3, and blue line corresponds to the moment t = 0.5. 




r r/(1+Kt) 

Fig. A.2. Snapshots of the evolution of the magnetic field components. Panel a) represents the de- 
pendence of the azimuthal component of the magnetic field on the radial coordinate for four time 
moments; Panel b) shows the dependence of the z-component of the magnetic field on the radial 
coordinate for four time moments; Panel c) and Panel d) illustrate the dependence of the modi- 
fied azimuthal and z-components of the magnetic field, respectively, on the self similar variable. 
Parameter values for this case are k — 2, Bo = 1 , Qout = 3£>in> Pin =0.1, p out = 1, ro = 1, and Oo = 1 . 
Black line corresponds to t — 0, red line represents the moment t = 0.5/3, the line shows the time 
moment 1/3, and blue line corresponds to the moment t - 0.5. 



